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Abstract 

We adress the question of how a c — c-state (a J/ip) can be dissociated by the 
strong color electric fields when moving through a color electric flux tube. The 
color electric flux tube and the dissociation of the heavy quarkonia state are both 
described within the Friedberg-Lee color dielectric model. We speculate on the 
importance of such an effect with respect to the observed J/'i/>-suppression in 
ultrarelativistic heavy ion collisions. 
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I. INTRODUCTION AND MOTIVATION 


Since the original work of Matsui and Satz who proposed J/V'-suppression as a clear signal 
for quark gluon plasma formation (the electric forces responsible for the binding are sufficiently 
screened by the plasma so that no boundstate should form) an intense work by several exper¬ 
imental groups (NA38-collaboration) has adressed this issue 0. Indeed a suppression in the 
corresponding dilepton signal compared to the Drell-Yan background has been seen over the last 
years for lighter projectiles like 0-l-Cu, 0-l-U and S-l-U. However, these observations may also 
be explained alternatively: (1) A part or even the total suppression can probably be explained 
by J/V’-absorption on the surrounding nucleons 0, i.e. J/ijj + N ^ A^ + D; (2) additional 
absorption by exothermic reactions like i.e. J/'ip + p^D + D might be attributed to ‘comovers’ 
(‘mesons’) being produced as secondaries Indeed, if ~ 7 mb is assumed, the reported 
suppression could be nicely reproduced by the absorption within the nuclear environment. In any 
case, the general consensus has then been that a highly (energy) dense intermediate reaction zone 
was needed to explain the observed suppression. 

Recently the data taken in Pb-|-Pb reactions at 160 AGeV by the NA50-collaboration [Q have 
given new excitement to the whole issue as a stronger absorption has been found than suggested 
by the models of absorption models on nucleons. It was immediately speculated @] that a large 
inner region of the reaction zone must absorb all J/'ifj (being ‘black’) in order to be compatible 
with observation. Though it was reported by Gavin and Vogt 0 and by Gassing and Ko ^ 
that the comover model also could explain the stronger suppression as relatively more mesons 
are produced for large Et events (or fluctuations) than for lighter systems, the latter calculations 
also show that this suppression takes place at energy densities that are much higher than those 
in which normal hadrons are expected to survive. 

In popular microscopic models, simulating the whole reaction of the collision, in the hrst few 
moments interaction strings are produced which subsequently decay (fragment) into secondaries 
(mesons, baryons, string-like hadronic resonances). These are best visualized as color electric 
flux tubes. Their creation is thought to happen by the exchange of a color octet gluon or by 
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colorless momentum transfer among target and projectile nucleons. (The color sources in the 
endcaps should be thought as a quark on the one side and an anti-quark on the other side, or as a 
quark and a di-quark, respectively.) These strings might also overlap (if their density is large) to 
form higher charged tubes, so-called color ropes . In any case, before a further hadronic-like or 
partonic-like state of matter does form, such a temporary build up and decay of strings is assumed 
to describe the very early collision phase. In fact one might call such an environment a precursor 
of quark matter. 


The J/il) (or better cc-state which later will form the J/V') as a rather heavy hadronic particle is 
also produced at the earliest state of the reaction in a hard collision among the nucleons. Thus it 
is natural to ask what happens if a (pre-) J/'^-state moves (or enters) into an environment of color 
electric strings. As the string carries a lot of internal energy (to produce the later secondaries) 
the quarkonia state might get absorbed and completely dissociated by the intense color electric 
field inside a single flux tube. This intuitive idea is the basic reason for our work. 


The dynamics of the flux tube we describe within the Friedberg-Lee model |1l0|,[II| by a semi- 
classical transport algorithm (section II). As a hrst step, within this model, we generate a bound 
cc-state with reasonable properties in order to mimic a. J/xjj (section III). We then adress how 
such a state gets dissociated when inserted into a chromoelectric flux tube and how the tube 
subsequently does split (section IV). We conclude with a summary and the implications of our 
investigation on the issue of JZ-^-suppression. 


II. DYNAMICAL REALIZATION OF THE FRIEDBERG-LEE MODEL 


Recently we have written down semiclassical transport equations for quark distributions IjT^pR 
starting from the phenomenological Lagrangian of Friedberg and Lee 

C = - mo - joff)'!' + -U((y)- , (1) 


where T denotes the quark helds, a is the color singlet scalar held representing the long range 
and nonabelian effects (of multi gluon exchange), and the last term contains the interaction of 
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the residual classical and abelian color fields Af,. All the nonabelian effects are assumed to be 
absorbed in the color dielectric function K{a) which is chosen such that k vanishes as a approaches 
its vacuum value a = outside the bag and k = 1 inside, 

k((T) = |1- —j^&(o-vac-0-). (2) 

^vac 

U{cf) = + B (3) 


is the self-interaction potential for the scalar cr-field containing cubic and quartic terms. Details of 
the numerical realization and the further approximations made are found in |]T3| , p!^ . One obtains 
two transport equations for the phase space distributions of quarks (/) and antiquarks (/) 


- m*{d^m*)d!^)f{x,p) = g^Pf,F>^‘'dPf{x,p) 

- m*{d^m*)d!^)f{x,p) = -g^p^F^^''dPj{x,p) , (4) 


which are numerically solved by applying the test-paticle method, a classical equation of motion 
for the scalar soliton field 


d,d^a + U'{a) + -^K'{a)F;,Fr + = 0 , 


( 5 ) 


which is solved using a staggered leapfrog algorithm |^, and an equation of motion for the 
confined colourelectric field D = kE 


V(kVAo) = -ja , 


( 6 ) 


which is determined by a two-dimensional finite element method (appropriate for cylinder sym¬ 


metrical configurations) developed by Mitchell |T5 


First we note that here we explicitly suppress the color magnetic fields, i.e. 2E‘^. (In 

a string like configuration with nearly constant radius Wilets and Puff have shown that due to 
the displacement current in the Maxwell equation magnetic fields are not produced |^.) / and 
/ describe the distribution of charged quarks and anticharged antiquarks (or diquarks) so that 
the colour charge density is expressed by 
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(7) 


Jo = 


(27r) 


£p{f{x,p) - f{x,p)) , 


whereas the scalar density ps = tl/tl/ entering into eq. (j^) reads 

P‘ = 7 ;^ f i‘p’^(f(x,p) + f(x,p)) . 


( 27 r) 


UJ 


( 8 ) 


(?7 = 4 accounts for spin and flavour degeneracy.) 


Within this approach one succeeds in simulating dynamically conhnement as (a) a single charged 
quark is forced back into the conhguration by the dielectric displacement D if it tries to escape 
11]; and (b) the dielectric displacement D = kE stays also within the conhguration as it has to 
vanish in the outside (and nonperturbative) vacuum where k —> 0 0. 

Static properties of soliton-like nucleons and mesons are succesfully reproduced by an appro¬ 
priate hxing of the scalar potential U{a) and the coupling constant ps- The vector coupling 
as = = 1-92 is hxed by assuming an effective string constant r = 1 GeV/fm when simulat¬ 

ing a colourelectric hux tube [1^. Here the volume averaged electric held {Ez) takes a value of 1.6 


fm ^ with Ez^max = 2.2 fm ^ along the z-axis. The radius of the string results as approximately 


0.8 fm. 


III. J/V'- AND D,D-STATES 


As a hrst step, within our model, we want to describe semiclassically a bound cc-state (syn¬ 
onymously a J/'0) and also a bound cg-state (synonymously a D). We start from a light meson 
state described within our model (using the parameters of the model ||I2|-|I4|) and then adiabati- 
cally increase the mass of the initial light quark(s) {rrii^q ~ 10 MeV) to a hnal heavy quark mass 
‘ruj,, = 1.56 GeV (taken from ref. |l|), i.e. 


rriqit) = mi^q{t = 0) -F (m/,c - mi^q{t = 0))t/Ts.o. for 0 < f < Ts.o. 

(9) 

mq{t) = TUf^c for t > Ts.o. ■ 

The change of the meson state with mass of both (or one) quarks is thus described fully dynami¬ 
cally. In hgs. 1 we see how the bag size shrinks from Vcms ~ 0.66 fm to Vcms ~ 0.25 fm for the case 
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of a cc-groundstate and to Vcms ~ 0.40 fm for the case of a cg-groundstate for a switching on time 
of T^.o. = 20 fm/c. We checked onr resnlts for varions switching on times Ts,o. = 20,30,50 fm/c 
and fonnd that the state described does not depend on the parameter r^.o.. The shrinking of the 
states is in accordance with the expectation of any baglike model. When the quark masses are 
increased to their hnal value, the pressure the more heavy quarks assert on the surface drops too, 
thus leading to the shrinking of the whole bag until it stabilizes at a smaller radius where pressure 
equilibrium is reached again. The radius of the heavy meson states (the ‘J/t/’’- and ‘D’-meson) 
are in the right range of what one would expect from other quarkonia model. In addition, £g. 2 
shows the spatial quark distribution and the soliton solution a{r) as a function of the radius for 
the cc-state at the beginning and the end of the evolution according to (||). 

When determining the overall mass of the states we find for rric = 1.56 GeV 

mcc = 4.10GeV (10) 

rricq = 2.25 GeV . 


These numbers have to be compared to (a) the mass of the rjc'. nrir^^ = 2.94 GeV and the J/i)-. 
^j/ip = 3.097 GeV and (2) to the mass of the H-meson: mo = 1-89 GeV and the D*-. mo* = 2.01 
GeV. The cc-state generated within the model (and within the parameters chosen to describe the 
light mesons) are thus about 1 GeV too heavy, the cg-state about 0.4 GeV, whereas the difference 
2mcq — nice ~ 0.4 GeV is in quite good agreement to 2mo — mj/^ = 0.68 GeV. This discrepancy 
one can relax by choosing a smaller bare mass of the charm quark. This mass is rigorously not 
known and only extracted by various charmonium models and might vary between 1 to 1.6 GeV 
p!8[] . Taking the a priori unknown charm quark mass as ~ 1.25 GeV (compare hgs. 1), the mass 
of the cc-state would lie slightly above the experimental value by about ~ 200 MeV, whereas the 
mass of the cg-state is lowered by about the same value, thus providing a much better description 
within our model. In the following we stay within our model states (0) (described by a bare 
mass of me = 1.56 GeV) as our conclusions will not depend on their masses and the results given 
in the next section will only weakly be affected by assuming a different mass for me within the 
afore mentioned range. 
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IV. J/^/;-DISSOCIATION BY A COLOR FLUX TUBE 


We now turn to the main issue of the present work: What happens to a cc-state entering the 
region spanned by a single chromoelectric flux tube? 

We prepare a chromodielectric flux tube (see also |^) by pulling the opposite charge distribu¬ 
tions of a (light) groundstate meson steadily apart. The quarks are propagated with a constant 
velocity up to a hnal spatial extent of about 8 fm (see e.g. £g. 3) within a time interval of about 
13 fm/c. At that given time to ~ 13 fm/c we keep the quark distributions at the endcaps fixed, 
thus providing the sources for the chromodielectric field. We then insert the ‘by hand’ 

into the central interior of the (by now) stationary flux tube, where the initial momentum and 
spatial distributions of the heavy quarks are chosen as the one obtained in the previous section. 
Initializing the cylindrically symmetric configuration in this way, we proceed by solving the full 
set of dynamical equations of motion, i.e. the two Vlasov equations (^) for the heavy quark and 
antiquark distribution, the equation (^ for the soliton field a and the one for the electric potential 


In £g. 3 we show the evolution of the total electric field energy as a function of time. As 
one notices, after inserting the bound state, the field energy decreases approximately linearly 
with time. What happens is that the heavy quark (described by its phase space distribution) 
is steadily pulled to one endcap of the string carrying opposite charge, whereas the antiquark 
carrying the opposite charge is correspondingly pulled into the other direction (see figs. 4). Being 
pulled apart, the two displaced charges screen the overall electric held in between. In other words, 
the electric held energy originally stored is transformed into the kinetic energy of the oppositely 
moving two heavy quarks. Being accelerated, after some smaller time interval At 1 fm/c, both 
heavy quarks have nearly reached the speed of light, so that the linear decrease of the overall 
electric held energy depicted in hg. 3 becomes obvious. 

In hgs. 4 also the evolution of the cr-held is depicted for various time stages. At the initial 
times one clearly sees how the soliton solution of the hux tube is distorted by introducing the 
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cc-state into the center of the string. At the timesteps up to about 18 fni/c the cr-£eld follows the 
quark distribution quite well, its shape reflects the shape of the latter distribution. At 20.8 fni/c 
the inner heavy charges have reached the hxed endpoints of the flux tube. At this point we stop 
- for numerical reasons - the further propagation of the quarks and concentrate on the further 
time development of the a-held shown in hgs. 5. In this further evolution the soliton held relaxes 
much more slowly to its vacuum value compared to the more or less instantanous screening of 
the electric held. It typically takes about 5 fm/c until a turns for the hrst time to its vacuum 


value JT^. This long time interval is basically a consequence of the nonlinear potential U{a) being 
hat around a ~ 0 and corresponds to the time needed for the soliton held to ‘move down’ the 
potential hill to its vacuum value at a = ay- Even when ay has been reached, though, the held 
still starts to oscillate around this value which can be seen in hgs. 5: In the second row, at t=23.2 
fm/c, the (T-held value is slightly above cr^ao whereas in the third row, at t=26.8 fm/c, its value is 
nearly in the perturbative region. The original held energy carried by the soliton held (the ‘Bag’ 
energy) being still present after the heavy quarks have already moved apart is initially more or 
less conserved. The late oscillations observed in the simulations are only damped by a transverse 
expansion of the held itself, which, however, is a rather slow process because of the large ‘glueball’ 
mass of the held. Thus it takes up to ~ 15 fm/c until the cr-held has finally relaxed and the 
static heavy D/D-meson states are formed. If one would consider a chiral 0(4)- extension for the 
soliton held one expects that these oscillations might accordingly transform into low-momentum 
(and nearly massless) pion modes and should accordingly be damped away more quickly. 


From our simulation we conclude that a localized cc-state immediately gets separated by the 
strong colourelectric held, and, in return, will finally be dissociated into H-meson (or A^-baryon) 
like configurations. One wonders if such a dissociation would persist also if initially the cc- 
state enters the string with some moderate transverse momentum. The hnite element method 
implemented in order to calculate the chromoelectric potential Ao{rt, z) can presently only handle 
cylindrically symmetric conhgurations (as e.g. a static string). Hence we are yet not able to 
simulate the full dynamics when a cc-state (as generated in the previous section) enters a static 
and elementary hux tube with some transverse momenta and hnite impact (though ths work is 
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in preparation |^). In the following we give an estimate that such a disintegration of a cc-state 
should still persist for even quite large transverse momenta. 


The time needed to pass the whole flux tube is roughly estimated as 


At 


Or 

r 


j/ip 




Vt 


\Jpt + 

\pt\ 


2 — 4fm/c 


( 11 ) 


for moderate to large (transverse) momenta in the range from 1 to 5 GeV/c. For an average held 
strength {\Ez\) of 1.61/fm^ the heavy quarks will suffer a gain in longitudinal momentum given 
approximately by (a^ = gl/Air ^ 2) 


\Apz\ ^ gy\Ez\At 


1.6 ^At 

tm 


3-6 GeV/c , 


( 12 ) 


whereas its internal rise in longitudinal kinetic energy of both quarks (by assuming that each 
heavy quark moves, on the average, with, let us say, half speed of light) 

AW ~ 2gy\Ez\As ~ gy\Ez\At ^3 — 6 GeV . (13) 


Such an increase in momenta or internal kinetic energy is certainly sufficient to break up the cc- 
state into two D-meson like states, as the energy needed is given by the difference in mass, i.e. of 
about 0.68 GeV. As the strings are thought to be energetic enough to produce all secondaries over 
the whole rapidity range in the heavy ion collision, its internal color electric held should be large 
enough to break up even more heavily bound systems like the J/'ifj. We are hence quite conhdent 
that our observed dissociation would also be found within a full 3-dimensional simulation of a 
cc-state interpenetrating a chromodielectric hux tube with a nonvanishing transverse momentum. 

In some popular charmonium models a rather strong ehective Goulomb potential between the c 
and c quark of the form We ~ —C/r , to be valid at short distances, gives a dominant contribution 
to the overall spectra and gives rise to a binding energy estimated to be 800 MeV for the J/ijj 
as a Is-state. When modeling this state within our semiclassical approach, the groundstate is 
only bound by the strong conhning scalar held. The dissociation described above corresponds to 
classical held ionization. The same ehects exist in a quantum mechanical treatment: A quantum 






mechanical system of two heavy quarks will dissociate at long times if a constant background 


held with a potential V ~ —g.uEz is switched on E). For our specihc numbers, along the z- 
direction, this leads to a potential of Id ~ (—1.7)zGeV/fm — (—2.1)2; GeV/fm due to the strong 
electric held. The modihcations are thus strong and nonperturhative. For 2 ; = 0.3 fm the ehective 
total potential Veff{z) ~ —Cjz — QvEz is already lowered by ~ 600 MeV compared to the ‘bare’ 
potential —Cjz. The existence of a (then highly modihed) resonant Is-state seems very unlikely. 
In addition one also can consider this quantum mechanical problem as a nonstationary situation 
when the J/V’ passes through the tube with some transverse velocity. When entering, the abrupt 
change in situation can not be considered as adiabatic, but better as diabatic: The original bound 
state has no time to adiabatically lower its energy by at least 0.6 GeV, but will immediately 
dissolve into the continuum. 


V. CONCLUSION AND OUTLOOK 


In this work we have adressed the question of how a cc-state (a T/V') behaves inside a chro¬ 
moelectric flux tube. The flux tube as well as the cc-groundstate are described within the model 


by Friedberg and Lee |^. Due to the strong color electric held inside the hux tube the heavy 
meson becomes dissociated rather immediately (on a timescale of 1 fm/c) ending up hnally in 
D- or Ac-like states. We could not simulate a fully 3-dimensional situation when the heavy meson 
enters the string with some hxed transverse momentum because our present numerical realization 
of the model works only for cylindrically symmetric conhgurations. However, we give an estimate 
that such a dissociation should also occur by held ionization for moderate and higher transverse 
momenta of the cc-state when passing through the tube within our semiclassical approach. 

We believe that this investigation adresses an intriguing question for understanding particular 
topics of QGD with respect to relativistic heavy ion collisions. If one believes in the succesful 
descriptions of microscopic approaches (like FRITIOF |^, RQMD |^, VENUS [^, HSD p^), 
a large region in space in the hrst few moments after the reaction is spanned by highly excited 
longitudinal strings. These strings will later materialize into all the secondaries and the total 
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energy carried by them, thus leading to a more or less complete distribution of the energy be¬ 
tween target and projectile rapidity. In this sense the flux tubes should be interpreted as highly 
energetic excitations of the QCD vacuum. One might consider such regions as a precursor of a 
QGP and deconhnement. Especially for the more heavy systems the effective region (or volume) 
of all the strings being produced within a short time interval of less then one fm/c gets so large 
that the strings become closely packed or already overlap. Although most often the strings are 
thought to fragment independently one might also consider the possibility of color rope formation 
as higher charged tubes. For a quantitative estimate we depict in fig. 6 the position of nucleonic 
strings (highly energetic strings excited by a target and projectile nucleon) all being produced 
in a cms-time interval of ~ 0.5 fm/c within a slab of the complete transverse area and 1 fm 
thickness along the longitudinal direction in the cms frame in a central S-l-U reaction at 200AGeV 
generated within the HSD approach . In total about 58 strings are produced in the transverse 
area where the sulphur ion faces the uranium nucleus. This area is about 40 fm^. If one assumes 
that the strings are homogeneously distributed (which they are not) this corresponds to a mean 
area of 0.7 fm^/string. Hence all strings are closely packed if the radius of a string is about 0.5 fm 
(in the Friedberg-Lee model it turns out to be nearly 0.8 fm). This simple reasoning illustrates 
that the whole reaction area in an ultrarelativistic heavy ion collisions might be completely hlled 
by strings in the first few moments. If such a large spatial region does exist, a large fraction of 
the produced J/V^’s have to pass it and thus can be affected by this highly excited environment. 
This will depend crucially on the (average) length of the produced strings before they hadronize, 
as this decides about the available ‘empty’ space for the (or pre-J/V'-states) to escape the 

initial reaction zone without entering any of the individual strings being built up. If they will 
get dissociated (in part or all of them) this would lead to a minor or stronger suppression of 
J/'ijj to he finally observed. Such a conclusion can only be further tested and strengthened if one 
incorporates these ideas within one of the present microscopic transport algorithms. 
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FIG. 2. The radial distribution of the soliton field and the quark density of the ‘heavy’ quark of the 
cc-state is depicted at the initial time (when the constituent quarks are light) and at a very late time 
(when the masses of the quarks have already incresed to their final value). 
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FIG. 3. The temporal evolution of the total electric field energy inside the string is shown. Up to 
times t ~ 13 fm/c the size of the string is linearly increased to a final length of ~ 8 fm and held constant 
afterwards. At that time t = 13 fm/c the cc-state is inserted into the central region of the string. A 
nearly linear decrease in time is observed. 
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FIG. 4. The quark distribution (left side) and the scalar field distribution (right side) at various times 
in the evolution (t = 12.8, 13.6, 16.2, 17.4 fm/c) are shown in a contourplot. 
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FIG. 6. The position of nucleonic strings is plotted, all being produced in a cms-time intervall of 
~ 0.5 fm/c within a slab of the complete transverse area and 1 fm thickness along the longitudinal 
direction in the cms frame in a central S+U reaction at 200AGeV generated within the HSD algorithm 
1^ . In addition, each string is accompanied by a circle of radius 0.5 fm. 
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